Analysis of Alternative mRNA Splicing in Vemurafenib-Resistant Melanoma Cells

Alternative splicing (AS) is one of the hallmarks of human cancer. One of the most common mechanisms of vemurafenib resistance in malignant melanoma is AS of BRAF, occurring in 15-30% of patients. The aim of our study was to investigate the transcriptome and AS D04landscape in the isogenic BRAF V600E cell line pair SK-MEL-239, where the vemurafenib-resistant derivative expresses a truncated BRAF transcript that lacks the RAS-binding domain. Transcriptome analysis showed differential expression of spliceosome components between the two cell lines. AS analysis, by four different tools, DEXSeq, rMATS, ASpli, and LeafCutter, has identified genes enriched for cell motility and melanin synthesis in vemurafenib-resistant cells. Overlapping predictions for all four tools have been experimentally validated. Our study expands the understanding of melanoma drug resistance from a new perspective and supports the need to investigate in detail the aberrant AS landscape in patients with malignant melanoma. Alternative mRNA splicing is common in cancers. In BRAF V600E-mutated malignant melanoma, a frequent mechanism of acquired resistance to BRAF inhibitors involves alternative splicing (AS) of BRAF. The resulting shortened BRAF protein constitutively dimerizes and conveys drug resistance. Here, we have analysed AS in SK-MEL-239 melanoma cells and a BRAF inhibitor (vemurafenib)-resistant derivative that expresses an AS, shortened BRAF V600E transcript. Transcriptome analysis showed differential expression of spliceosome components between the two cell lines. As there is no consensus approach to analysing AS events, we used and compared four common AS softwares based on different principles, DEXSeq, rMATS, ASpli, and LeafCutter. Two of them correctly identified the BRAF V600E AS in the vemurafenib-resistant cells. Only 12 AS events were identified by all four softwares. Testing the AS predictions experimentally showed that these overlapping predictions are highly accurate. Interestingly, they identified AS caused alterations in the expression of melanin synthesis and cell migration genes in the vemurafenib-resistant cells. This analysis shows that combining different AS analysis approaches produces reliable results and meaningful, biologically testable hypotheses.Peer reviewe

Analysis of Alternative mRNA Splicing in Vemurafenib-Resistant Melanoma Cells

Alternative splicing (AS) is one of the hallmarks of human cancer. One of the most common mechanisms of vemurafenib resistance in malignant melanoma is AS of BRAF, occurring in 15-30% of patients. The aim of our study was to investigate the transcriptome and AS D04landscape in the isogenic BRAF V600E cell line pair SK-MEL-239, where the vemurafenib-resistant derivative expresses a truncated BRAF transcript that lacks the RAS-binding domain. Transcriptome analysis showed differential expression of spliceosome components between the two cell lines. AS analysis, by four different tools, DEXSeq, rMATS, ASpli, and LeafCutter, has identified genes enriched for cell motility and melanin synthesis in vemurafenib-resistant cells. Overlapping predictions for all four tools have been experimentally validated. Our study expands the understanding of melanoma drug resistance from a new perspective and supports the need to investigate in detail the aberrant AS landscape in patients with malignant melanoma. Alternative mRNA splicing is common in cancers. In BRAF V600E-mutated malignant melanoma, a frequent mechanism of acquired resistance to BRAF inhibitors involves alternative splicing (AS) of BRAF. The resulting shortened BRAF protein constitutively dimerizes and conveys drug resistance. Here, we have analysed AS in SK-MEL-239 melanoma cells and a BRAF inhibitor (vemurafenib)-resistant derivative that expresses an AS, shortened BRAF V600E transcript. Transcriptome analysis showed differential expression of spliceosome components between the two cell lines. As there is no consensus approach to analysing AS events, we used and compared four common AS softwares based on different principles, DEXSeq, rMATS, ASpli, and LeafCutter. Two of them correctly identified the BRAF V600E AS in the vemurafenib-resistant cells. Only 12 AS events were identified by all four softwares. Testing the AS predictions experimentally showed that these overlapping predictions are highly accurate. Interestingly, they identified AS caused alterations in the expression of melanin synthesis and cell migration genes in the vemurafenib-resistant cells. This analysis shows that combining different AS analysis approaches produces reliable results and meaningful, biologically testable hypotheses.Peer reviewe

Signalling ballet in space and time.

Although we have amassed extensive catalogues of signalling network components, our understanding of the spatiotemporal control of emergent network structures has lagged behind. Dynamic behaviour is starting to be explored throughout the genome, but analysis of spatial behaviours is still confined to individual proteins. The challenge is to reveal how cells integrate temporal and spatial information to determine specific biological functions. Key findings are the discovery of molecular signalling machines such as Ras nanoclusters, spatial activity gradients and flexible network circuitries that involve transcriptional feedback. They reveal design principles of spatiotemporal organization that are crucial for network function and cell fate decisions

Dynamics of receptor and protein transducer homodimerisation

<p>Abstract</p> <p>Background</p> <p>Signalling pathways are complex systems in which not only simple monomeric molecules interact, but also more complex structures that include constitutive or induced protein assemblies. In particular, the hetero-and homo-dimerisation of proteins is a commonly encountered motif in signalling pathways. Several authors have suggested in recent times that dimerisation relates to a series of physical and biological outcomes used by the cell in the regulation of signal transduction.</p> <p>Results</p> <p>In this paper we investigate the role of homodimerisation in receptor-protein transducer interactions. Towards this end, mathematical modelling is used to analyse the features of such kind of interactions and to predict the behaviour of the system under different experimental conditions. A kinetic model in which the interaction between homodimers provokes a dual mechanism of activation (single and double protein transducer activation at the same time) is proposed. In addition, we analyse under which conditions the use of a power-law representation for the system is useful. Furthermore, we investigate the dynamical consequences of this dual mechanism and compare the performance of the system in different simulated experimental conditions.</p> <p>Conclusion</p> <p>The analysis of our mathematical model suggests that in receptor-protein interacting systems with dual mechanism there may be a shift between double and single activation in a way that intense double protein transducer activation could initiate and dominate the signal in the short term (getting a fast intense signal), while single protein activation could control the system in the medium and long term (when input signal is weaker and decreases slowly). Our investigation suggests that homodimerisation and oligomerisation are mechanisms used to enhance and regulate the dynamic properties of the initial steps in signalling pathways.</p

The secret life of kinases: functions beyond catalysis

Protein phosphorylation participates in the regulation of all fundamental biological processes, and protein kinases have been intensively studied. However, while the focus was on catalytic activities, accumulating evidence suggests that non-catalytic properties of protein kinases are essential, and in some cases even sufficient for their functions. These non-catalytic functions include the scaffolding of protein complexes, the competition for protein interactions, allosteric effects on other enzymes, subcellular targeting, and DNA binding. This rich repertoire often is used to coordinate phosphorylation events and enhance the specificity of substrate phosphorylation, but also can adopt functions that do not rely on kinase activity. Here, we discuss such kinase independent functions of protein and lipid kinases focussing on kinases that play a role in the regulation of cell proliferation, differentiation, apoptosis, and motility

The Renilla luciferase gene as a reference gene for normalization of gene expression in transiently transfected cells

Background: The importance of appropriate normalization controls in quantitative real-time polymerase chain reaction (qPCR) experiments has become more apparent as the number of biological studies using this methodology has increased. In developing a system to study gene expression from transiently transfected plasmids, it became clear that normalization using chromosomally encoded genes is not ideal, at it does not take into account the transfection efficiency and the significantly lower expression levels of the plasmids. We have developed and validated a normalization method for qPCR using a co-transfected plasmid. Results: The best chromosomal gene for normalization in the presence of the transcriptional activators used in this study, cadmium, dexamethasone, forskolin and phorbol-12-myristate 13-acetate was first identified. qPCR data was analyzed using geNorm, Normfinder and BestKeeper. Each software application was found to rank the normalization controls differently with no clear correlation. Including a co-transfected plasmid encoding the Renilla luciferase gene (Rluc) in this analysis showed that its calculated stability was not as good as the optimised chromosomal genes, most likely as a result of the lower expression levels and transfection variability. Finally, we validated these analyses by testing two chromosomal genes (B2M and ActB) and a co-transfected gene (Rluc) under biological conditions. When analyzing co-transfected plasmids, Rluc normalization gave the smallest errors compared to the chromosomal reference genes. Conclusions: Our data demonstrates that transfected Rluc is the most appropriate normalization reference gene for transient transfection qPCR analysis; it significantly reduces the standard deviation within biological experiments as it takes into account the transfection efficiencies and has easily controllable expression levels. This improves reproducibility, data validity and most importantly, enables accurate interpretation of qPCR data

When kinases meet mathematics: the systems biology of MAPK signalling

The mitogen activated protein kinase/extracellular signal regulated kinase pathway regulates fundamental cellular function such as cell proliferation, survival, differentiation and motility, raising the question how these diverse functions are specified and coordinated. They are encoded through the activation kinetics of the pathway, a multitude of feedback loops, scaffold proteins, subcellular compartmentalisation, and crosstalk with other pathways. These regulatory motifs alone or in combination can generate a multitude of complex behaviour. Systems biology tries to decode this complexity through mathematical modelling and prediction in order to gain a deeper insight into the inner works of signalling networks

Basic fibroblast growth factor modifies the hypoxic response of human bone marrow stromal cells by ERK-mediated enhancement of HIF-1α activity

Human bone marrow stromal cells (hBMSCs, also known as bone marrow-derived mesenchymal stem cells) are promising tools for the cellular therapy of human pathologies related to various forms of hypoxia. Although the current concepts of their clinical use include the expansion of hBMSC in standard cell culture conditions, the effect of the mitogen-driven ex vivo expansion on the adaptation to the hypoxic environment is unknown. Here, we provide data that the basic fibroblast growth factor (FGF2) enhances the induction of a wide range of hypoxia-related adaptive genes in hypoxic hBMSCs. We identified that the FGF2 signal is transmitted by the ERK pathway similar to that of hypoxia that also utilises the distal elements of the same signalling machinery including the extracellular signal-regulated kinase 1/2 (ERK1/2) and mitogen-activated protein kinase kinases (MEK1/2) in hBMSCs. We found that the simultaneous activation of ERK1/2 by FGF2 and hypoxia transforms the activation dynamics from oscillatory into sustained one. Activated ERKs co-localise with stabilised hypoxia inducible factor-1α (HIF-1α) followed by the reduction of its nuclear mobility as well as increased DNA binding capacity leading to the up-regulation of hypoxia-adaptive genes. Our findings indicate that the status of the ERK pathway has significant impacts on the molecular adaptation of hBMSCs to the hypoxic milieu